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I. INTRODUCTION 

A central problem of strong interaction theory is to derive the nucleon-nucleon interaction 
from more fundamental principles. It is well established that quantum chromodynamics 
(QCD) is the correct theory of strong interactions. The phenomenology of the nucleon- 
nucleon interaction suggests a description in terms of nucleons and the exchange of mesons, 
while QCD is a non-abelian gauge theory of quarks of three colors interacting by exchange of 
gluons. These two pictures, one of quarks and gluons and the other of nucleons exchanging 
mesons are difficult to reconcile. The notions of asymptotic freedom and its counterpart 
confinement suggest that quarks and gluons are the appropriate degrees of freedom for short 
distances and high energy, but are not for the large distances and relatively low energies 
of the nucleon interaction, where the quark picture becomes strongly non-perturbative. 
The problem is to make a theory, starting from QCD, that obtains the nucleons and their 
meson exchange. If such a theory agrees with the vast body of nucleon-nucleon interaction 
phenomenology, it will have, in some sense, "derived" traditional nuclear physics from QCD. 
This paper reports the details of an effort based on the Skyrme model [^^ to make that 
connection. A brief summary of these results have been presented elsewhere 0. 

The Skyrme model describes a nonlinear, classical field theory of interacting pions. This 
seems far from QCD, but in fact the Skyrme Lagrangian has been shown to be an effective 
theory of QCD in the large distance (long wave length) non-perturbative regime in the limit 
of a large number of colors (A^c)- There is no unique path from such a classical theory back 
to quantum chromodynamics with three colors, but attempts to construct alternate classical 
pictures have shown that there is surprisingly little leeway 0, and work to extract nucleon 
phenomenology from the Skyrme Lagrangian has been relatively successful with results that 
are fairly insensitive to variations in details of the theory. Hence the Skyrme Lagrangian 
provides a robust, QCD-based starting point for studies of nucleon and nucleon-nucleon 
properties in the long wavelength, non-perturbative regime. 

The Skyrme Lagrangian possesses a topologically conserved charge, B, that Skyrme |I| 
identified with baryon number. The B = 1 solutions of the Skyrme Lagrangian yield an 
appealing picture of the nucleon and its excited states, but since the theory is a large A'c 
approximation, corrections for real nucleon observables are expected to be of order 1/A^c 
or ~ 30%. An interesting and important problem area is the calculation of these 1/Nc 
corrections. These are first quantum corrections, and their calculation faces the standard 
uncertainties found in going from any classical theory, particularly a non-linear field theory, 
to a quantum theory. 

The relative success of the Skyrme picture for the B = 1 system invites its extension to 
B = 2 and in particular to the problem of the nucleon-nucleon interaction. The phenomenol- 
ogy of the A^A^ interaction suggests that, with the possible exception of the repulsive core, 
the interaction is well described in terms of multiple pion exchange. Thus the Skyrme 
model which expresses the physics in terms of interacting pions, seems to be a good starting 
point. Skyrme himself [|I| noted that for very large separation, the A^A^ interaction reduces 
to one pion exchange in his model. The key question then is whether the model can also 
give the central mid-range attraction that is responsible for nuclear binding. There are 
three problems that must be overcome in addressing this question. The first is the solution 
of the classical B = 2 field equations to obtain and describe B = 2 configurations. The 



second is the problem of obtaining nucleon-nucleon interactions from Skyrmion-Skyrmion 
interactions. The third is the introduction of 1/Nq corrections. As we shall see the last 
two are inexorably mixed. The need for careful 1/Nq corrections is much greater in the 
B = 2 problem than for B = 1, because the only way presently known to calculate the 
interaction energy is to calculate the energy of the total B = 2 system and then subtract 
the energy of two free B = 1 Skyrmions. This is a difficult, delicate, and dangerous process 
since the interaction energy is a very small fraction of the total. In this paper we present a 
procedure for carrying out the three steps. Our procedure yields a A^A^ interaction closely 
resembling standard phenomenology. The interaction has a one pion exchange tail at long 
distance, a short range repulsion, a tensor force and, best of all, a strong, mid-range central 
attraction-all in substantive agreement with phenomenology. 

The problem of solving or even of describing the B = 2 interacting Skyrmion system 
is complex. There are only two known stable (time-independent) B = 2 configurations. 
One is the trivial one of two infinitely separated B = 1 Skyrmions, and the other is a 
compact bound state of toroidal baryon density in which the identity of the two individual 
Skyrmions is lost. To obtain an interaction energy one must define a set of intermediate, 
unstable, static Skyrmion configurations at finite separation. Each configuration defines a 
point on an adiabatic, six dimensional, collective manifold, called the unstable manifold 
by Manton [^. Each point on the manifold is stabilized by imposing a constraint. The 
dimensionality of this manifold makes it virtually impossible to calculate the interaction 
energy at every point. Fortunately much of the interaction energy can be described by 
studying only a very small part of that manifold. But on that sub-manifold the description 
must be accurate. As we emphasized above, the interaction energy is calculated as the 
difference between large numbers. Hence seemingly small inaccuracies in configurations can 
lead to large discrepancies in interaction energy. It is this problem that has plagued early 
attempts to study the A^A^ interaction in the Skyrme model with the product ansatz. This 
ansatz, also first given by Skyrme |I]], describes the B = 2 system as the product of two 
B = 1 systems. The configurations so obtained are not solutions of the constrained non- 
linear field equations, and hence the energy they yield is above the true energy. For very 
large separation these differences are not serious, but for intermediate and short distances 
they are fatal. For example the product ansatz cannot describe the toroidal lowest B = 2 
state. Early efforts [^,0] to study the A^A^ system with the product ansatz developed much 
of the vocabulary and sense of the problem, but in the end failed to give any mid-range 
attraction. Surprisingly enough one finds a repulsive core in this ansatz, but it has been 
shown that this does not have the correct g-parity [§. A far better approach to the B = 2 
system is direct numerical solution of the constrained, non-linear, field equations. Recently 
Walhout and Wambach (WW) have produced just such solutions. They form the starting 
point of our work. 

How does one go from the Skyrmion-Skyrmion interaction energy calculated on the 
interesting sub-manifold of the unstable manifold to the A^A^ interaction? A commonly 
used prescription is to project the Skyrmions to nucleons at each separation. This method 
has been used extensively in the early product ansatz work. Using the same approach 
WW find weak central attraction. The nucleons only requirement is too strong, however. 
Phenomenologically all that is required is that the particles be nucleons asymptotically. 
When they approach and interact they can deform and be excited, for instance to the A 



state, as the dynamics dictates. This state mixing must be very important, since in the 
solutions for the interacting Skyrmions found by WW, the individual Skyrmions deform 
strongly as they interact, finally loosing their individuality completely in the toroidal state. 
If the time scales (or equivalently energy scales) for the radial motion and the excitation of 
Deltas separate, the contribution of this distortion to the energy can be calculated in the 
Born-Oppenheimer approximation. Finite Nq effects are critical in making the time scales 
appropriately distinct. It is this combination of finite A'"c effects and Born-Oppenheimer that 



is crucial to obtaining the phenomenological potential. Each idea has been used before |]T0 
but it is the combination of exact solution of the classical field equations, state mixing using 
the Born-Oppenheimer approximation and finite A'"c effects that together makes physical 
sense and gives the phenomenologically correct potential as we shall see below. We shall 
first treat the effects of the distortions perturbatively, where it is clear that the distortion 
makes an attractive contribution. Since the contribution of the distortions becomes large at 
small distances, we turn to diagonalization of the Born-Oppenheimer Hamiltonian to check 
the perturbative results. Once we make the Born-Oppenheimer separation we can calculate 
explicit corrections to the adiabatic picture. These take the from of transition potentials 
between the different adiabatic channels, as well as some velocity dependent potentials 
similar to Berry's phase terms. 

In Section II we set out some key notions. We review the Skyrme Lagrangian, the 
adiabatic manifold for interacting B = 2 Skyrmion systems, review the algebraic method for 
projecting baryon states from Skyrmions and for introducing finite Nq effects, and show how 
these algebraic methods permit a description of the Skyrmion-Skyrmion interaction on the 
unstable manifold. This section brings together the known tools we need for the problem. 
Section III presents our derivation of the NN interaction from the Skyrmion interaction for 
finite Nq. It gives the adiabatic interaction derived first in perturbation theory and then 
using the Born-Oppenheimer approximation. The formalism for obtaining the transition 
potentials is also given here. Our results are presented in Section IV. As we have discussed 
above, the potentials are close to static A^A^ potentials obtained from the data. Section V 
gives a brief summary and conclusions. Some of the algebraic details of the perturbation 
treatment are presented in the appendix. 

II. SOME KEY NOTIONS 

A. The Skyrme model 

The Skyrme model is a non-linear field theory that can be realized in terms of an SU(2)- 
valued matrix field U, with Lagrangian density 

C = &TT[d,U{x)d^U^{x)] + ^TT[U^d,U, U^dMf. (1) 

The model is covariant, as well as invariant under global SU(2)-rotations that are identified 
with the isospin symmetry. As was discovered by Skyrme the model has a topologically 
conserved quantum number, which is identified as the baryon number B. The U field is 
interpreted as a combination of a scalar a field and an isovector pion field, U = a + if ■ -k. 
The a field is not an independt physical field due to the unitarity constraint on U . 



The standard time-independent solution to the classical field equations for B = 1 is the 
defensive hedgehog, where the pion field points radially outward, 

f/i (r) = exp(ir- r/(r)). (2) 

The baryon number of this state is given hj B = (/(O) — /(oo))/7r = 1. This solution breaks 
translational invariance, as well as the 0(4) spin-isospin symmetry. If we perform a global 
SU(2) isorotation on the state, 

U^ir\A) = A^U^ir)A, (3) 

we obtain a state of the same energy. One can easily show that an isorotation A has the 
same effect on U as performing an ordinary rotation i?[y4~^] on f, so that the grandspin K, 
the sum of spin and isospin, is a good quantum number, with eigenvalue 0, 

K = S + I. (4) 

In the B = 2 system we shall frequently use the product ansatz. This ansatz makes 
use of the fact that the product of two B = 1 solutions has baryon number two. The most 
general ansatz we can therefore construct from two hedgehogs consists of the product of two 
rotated hedgehogs, where the centers have been translated, 

U2ir\RAB) = A^Ui{f- R/2)AB^Ui{f+ R/2)B 
= U2{r\RCD) 
= D^C^'^^Ui{r- R/2)CU,{f+ R/2)C'/^^D. (5) 

In the last of (^ we have introduced the matrix D, that describes the rigid isorotation of 
the whole system, as well as a relative isorotation C. When R is very large changing C or 
D does not change the energy of the solution. For smaller R, D still generates a zero-mode 
(corresponding to broken isospin symmetry), but the energy will depend on C. Again, the 
energy is also invariant under spatial rotation, due to the conservation of angular momentum 
J = L + S. 



B. Determination of manifold 

After having constructed all the constrained B = 2 solutions, we construct a manifold 
that contains all these solutions. This method, which is well know in nuclear structure 



physics as the adiabatic large amplitude collective motion |11], seems first to have been 
introduced in the theory of interacting Skyrmions by Manton |^, under the name unstable 
manifold. His discussion is based on the deep geometric results by Atiyah and Hitchin for 



the interaction of BPS monopoles [|T2| , which make an identification of the structure of the 
manifold very clear but do not give many tools for constructing it. 

One of the ways to solve for the adiabatic collective manifold would be to solve a local 
fiuctuation equation. This is a very complex problem for a field theory, but even when 
we reduce the complication by selecting only a finite number of degrees of freedom as in 
the Atiyah-Manton ansatz |ir3| , p!4| this remains a formidable (but interesting) task. The 



minimal set of degrees of freedom can be identified without performing any calculation, 
however. There are six trivial ones, consisting of the three center of mass coordinates, as 
well as global rotations of the whole system (in both space and isospace), which we shall 
ignore in the present discussion. For large separations we find that an additional six degrees 
of freedom are necessary to describe the product ansatz (|^): the relative separation R (three 
coordinates) and three more coordinates describing the relative orientation (as discussed 
in the previous subsection usually represented by the U{2) matrix C = c^I + ic- t, where 
c1 + (^ = 1). The minimal degrees of freedom are just the continuously deformed set of 
six parameters as we take R to smaller values. The meaning of the word "adiabatic" in 
adiabatic manifold is that the motion in the manifold consists of the slowest modes. This 
only makes sense if at any point all the unstable directions lie within the collective space, 
so that no lower energy mode lies outside the manifold. An analysis of the unstable B = 2 
hedgehog by Manton f^ shows that it is plausible that this state has six unstable modes, 
again corresponding to relative separation (three modes) and a relative rotation (another 
three). Since this hedgehog is supposed to be part of the collective surface, we again find a 
minimum of six collective modes. 

Once one determines the collective surface there exist standard techniques |[ll[] to de- 



termine the collective Hamiltonian, with the restriction that it is quadratic in momenta. 
Schematically it takes the form 

H = K{R, R, C,pR, L,,pcj) + \liM,^Ij + \J^N^JJ, + V{R, R, C) (6) 

Here pr is the momentum associated with R, L the orbital angular momentum and pc is 
shorthand for the momentum associated with C. The symbols Jj and J, denote the intrinsic 
components of the (conserved) isospin and angular momentum, respectively. The kinetic 
energy is closely related to the geometric structure of the collective surface, and can be 
very complicated. One example that has been studied extensively is the case of the BPS 



monopole, where one has only kinetic and no potential energy. It has been proven |T^ 
that the underlying geometry is reasonably complicated and allows for non-trivial scattering 
of monopoles. In the current work we concentrate on the potential energy, however. We 
assume that the kinetic terms (including the contributions of / and J) take the form 

K^^fn^j^iCUcqi (7) 

where C2 denotes the 0(4) Casimir operator ^{Sf + If). The only correction to this that 
has been studied is the R and C dependence of M [^, but we ignore that in this paper. It 
will be interesting to return to this problem in the context of the NN interaction. 

To calculate the manifold beyond the product ansatz we use the fact that for special 
choices of C, due to symmetry, there exist three stable situations WM - the stable attractive 



channel, the marginally stable hedgehog-hedgehog and repulsive channels - where a path 
of the collective surface is traced out for fixed C. In each of these channels (labeled by 
i for convenience) one calculates the classical potential energy V^. The identification of 
the symbols R in each of the three Hamiltonians as one and the same radial coordinate 
is a bit precarious. There is no apparent reason why we cannot separately redefine the 
radial coordinate in each channel (say, by a point transformation). Here the topology comes 



into play however. Since each channel corresponds to a definite choice of C, the three point 
transformations can be embedded in a single transformation that depends on both R and C, 
so that we are just making a change of coordinate transformation on the manifold, by which 
we do not gain much. In order to sensibly identify the i?'s in the same channel we should 
use geometric notions. First note that each of the three channels i seems to correspond to a 
geodesic line on the collective surface, due to the extra symmetry of the pion field. In order 
to be able to make a meaningful identification of points on the three lines as having the same 
R the three points should lie on a single geodesic surface on the manifold, that crosses the 
three geodesic lines i orthogonally. This is of course extremely difficult to implement, and 
we do not know of a good way of doing this. In most discussions of the Skyrmion-Skyrmion 
interaction, including the present one, these problems are ignored, and a special definition 
of R is made that is then used to identify the different channels. 

C. The algebraic model 

Having determined the functions VJ we need to understand their relation to the nucleon- 
nucleon potential. The usual approach is to return to the product ansatz, where one makes 
very strong assumptions about the form of the two-Skyrmion solution. Again, this form is 
exact in the large-i? limit. The advantage of this form is that it guides us into a partic- 
ular expansion of V in terms of functions c ■ R, C4 (the quaternion c specifies the relative 
orientation). This expansion is closely related to an algebraic model of the Skyrmion ||15|| , 



which in the large- iV limit behaves in exactly this manner. There is no a priori reason why 
this form should be restricted to the product ansatz, however. As we have argued above 
the coordinates of the collective surface are in one-to-one correspondence with those of the 
product ansatz, so that the current problem can again be modeled by a similar form |1^ . 



Let us now discuss the algebraic model [ll8| , p!5| that we use. Originally it was introduced 



to give an algebraic alternative to the rotational zero-mode quantization of the Skyrmion. 
Since the Skyrmion breaks SO{A) ~ SUi{2) x SUs{2) symmetry, but is invariant under a 
5*0 (3) subgroup, we find that we can describe each of the equivalent configurations by a 
point of on the sphere ^3 ~ S0{A)/S0{3) ^. 

We can now introduce a bosonic model with dynamical symmetry SU{A), which in the 
limit of large boson number is also described by a point on the sphere 5*3. Without discussing 
the hamiltonian in any detail, we can draw some interesting conclusions. For any A^ (the 
number of bosons) we find that, as in the Skyrme model, all states have I = S. The 
spectrum of the algebraic model is cut off at / = N/2, however. Since we believe that in 
the real world, where A'c = 3, we have a cut-off at isospin 3/2 for the delta, we impose the 
condition that the number of bosons A^ is equal to the number of colors A'c. If we make 
this identification one can immediately show that for A^ = A^c = 3 the algebraic model is 
identical to the non-relativistic quark model without radial excitations. It has been shown 
that this model can be used to extract finite A'c corrections to observables. Of course we 
quantize only a few coordinates and we can only calculate a very limited set of finite A'c 
corrections, corresponding to the zero-mode quantization. One may hope, however, that 
through a judicious choice of collective coordinates we have found the most important, and 
maybe even the largest, contributions. 

In this paper we are not interested in the single Skyrmion case, but rather in the the case 



of two Skyrmions. This has been discussed in quite some detail for the case of the product 



ansatz in Ref. ^^. The algebraic model used consists of two sets of m(4) algebras, one for 
each Skyrmion, as well as a radial coordinate R. The interaction of the two Skyrmions can 
be expanded in terms of three operators, the identity /, and the operators W and Z: 

Z = T^;r^^[?>R,R,-5,,]/Nl (8) 

Here a and (3 label two different sets of bosons, used to realize the ^(4) algebras, and T is 
a one-body operator with spin and isospin 1. The semi-classical (large-A'^c) limit of these 
operators can be given in terms of R and C = C4I + c ■ r as |T5| : 

Wci = 3cl - c^, 
Zci = 6(c-i?)2-2c^. (9) 

If we expand the Skyrmion- Skyrmion interaction we find 

v{R) = vi{R) + V2{R)Wci + V3{R)Z^i + higher order. (10) 



Fortunately studies of the B = 2 system with the product ansatz [pi|7|p!5| , the Atiyah-Manton 
ansatz [0,0 and with direct numerical integration of the equations have shown that 
much information about the B = 2 energy as a function of separation R is contained in the 
particular choice of three paths discussed earlier, corresponding to three definite choices of 
the relative orientation, or equivalently to a special symmetry of the pion field. This allows 
us to calculate an approximate form for the interaction energy of two Skyrmions which 
consist of the first three term of (|1^) (the form used in Ref. 0] is fully equivalent), which 
correspond to an expansion up to first order only in terms of the two scalar operators W 
and Z. One of the important consequences of the explicit realization of these operators is 
that now algebraic methods [15,18| can be used to obtain their explicit finite A^c realization. 



In this way we can include some (arguably the most important) finite A^c corrections in 
the calculation. This includes the well-known amplification of the nucleon-nucleon matrix 
element of T by the factor {Nq + 2)/Nc. 

III. THE ADIABATIC INTERACTION 

The problem we address here is how to obtain a nucleon-nucleon interaction which in- 
cludes some finite A'c effects, from ([lOD . To date most workers have simply sandwiched 
(p!0|) between states containing two nucleons for each value of R and called that the interac- 
tion. This is definitely not the nucleon-nucleon interaction determined phenomenologically 
from scattering phaseshifts, where we can only require that the state contains two nucle- 
ons at large separation. For shorter distances, when the hadrons are interacting they can 
be (virtually) what ever the dynamics requires, for example A's. We know |]T^ that the 
baryon density is strongly deformed for moderate R - a. deformation that finally results in 
the toroidal shape with its 70 MeV of binding. That configuration corresponds to a great 
deal of mixing between all states of the individual Skyrmions. This admixture will lead to 
additional attraction at intermediate distances. The purpose of this paper is to calculate 
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the effect of this mixing on the N-N interaction for A'^c = 3 starting with ([T0|). The need 
for state mixing has been reahzed before [^ in the context of the product ansatz, but that 
ansatz is a poor starting point. Furthermore, as we state below, state mixing only makes 
sense for finite Nc- At the same time there are algebraic finite A^c corrections to matrix 



elements that are difficult to include in the formalism of ||2T|. For these reasons, these early 
attempts to introduce state mixing did not find sufficient central attraction. We should 
remark here that a very different formalism, based on a study of fluctuations around the 



product ansatz |^ has been used to show that some attraction exists in two-pion range. 
Unfortunately the expansion underlying this approach fails in the region that is most crucial 
to central attraction; 1 f m < i? < 2 fm. 

For large A'^c, one can distinguish two energy scales or reciprocally two time scales. The 
slower time scale is associated with the motion in the collective manifold, i.e., R and the 
orientation. The other time scale corresponds to the almost instantaneous response of the 
pion field to changes in R and the relative isospin orientation. For large Nq we cannot 
separate the time scales for the two sets of adiabatic modes, as can be seen in the highly 
correlated doughnut. For A'"c equal to three the situation changes. The R motion is typically 
much slower than the rotational motion which leads to the separation of the nucleon and A 
states. We thus have three energy scales, that of the pion field, of the A^ — A separation and 
of the R motion. We can now calculate a Born-Oppenheimer potential for the i?-motion, 
which constitutes the slowest degree of freedom. 

A. Pert ur bat ive approach 

All of the effects of quick response of the pion field at each R are already in Eq. (plQl). It 
is the effect of the rotational states we want to include. As a guide, we begin by studying 
it perturbatively. Let us call the full potential between two nucleons, including the effects 
of the rotational excited states, V. In terms of the Skyrmion interaction v of (|r^) and to 
second order, V is given by 

V(«)HA'A'K«)|A'A') + E'«»»^. (") 

(Here Ejy^ is the two-nucleon energy and Es is the energy of the relevant excited state.) 
The first term on the right is the direct nucleon-nucleon projection of v and is the term 
that has appeared in the literature. The second term is the correction due to rotational or 
excited states. It is clear from the energy denominator that the second term is attractive. 
For A'c = 3 the states \s) that can enter are |A^A), |AA^) and |AA). Except for a centrifugal 
term we shall add later, the excitation energy is assumed to be independent of R, and thus 
can be expressed in terms of the N-A energy difference, 300 MeV. It is here that we are 
neglecting the effect of the position dependence of the inertial parameters. 

The lowest order potential has been calculated by several groups using several levels of 
approximation to the Skyrme model [0,mH|. The form found (or assumed), however, can 



always be cast in the form of the algebraic model given by |T5|. This allows us to study both 
the large- A"c limit, as well as to include finite A"c effects explicitly in a systematic way. (The 
reader should note that the results in |^ include corrections for the tensor and spin-spin 
interactions.) 
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The leading term in this expansion is given by the form (cf. ([T0|)) 

v{R) = v^{r) + V2{R)W + v^{R)Z{R). (12) 

The algebraic operators W and Z have very simple expectation value for nucleons, (see |p!5[ ) 

(A^A^I W \NN) = \Pl {NN\ o^-a^T^- r^ \NN) , 
(A^A^I Z \NN) = iPJj {NN\ (3(T^- Ra^-R-a^- a'^^- r^ |A^A^) , (13) 

i.e., apart from a finite Nq correction factor P^ {Pisf = 1 + 2/Nc) they represent the spin- 
spin and tensor interactions. If we use this form we find that the lowest order interaction in 
the nucleon-nucleon channel is given by 

If we wish to evaluate the perturbation correction ( [TTl) we must evaluate the two sets 
of matrix elements {NN\ v{R) \NA) {NA\ v{R) \NN) and {NN\ v{R) |AA) (AA| v{R) \NN) 
separately. Details are discussed in appendix A, but the final result for the first order 
correction to the NN interaction is 



^PT = -^miPo + i¥rPN + 4fQ"N)Pnivl + 2vl) + 
(^'■^')h^sQNPo + HPn - i-,Ql)Pnivl - vl] 



{3a'- Ra'-R-a'- a')[-±QlP^^ + (ifp^ _ ^^Ql)P[]{vl - v^v,)}. (15) 

Here Qn is an additional finite A^c correction factor 



Qn = ,J{l-l/Nc){l + 5/Nc), (16) 

6 is the N-A energy difference, and Pf is a projection operator onto isospin T. 

B. Diagonalization of the interaction 

The perturbative expansion discussed in the previous subsection just gives a rough es- 
timate of the effect of state mixing in the Born-Oppenheimer approximation. To solve for 
the Born-Oppenheimer Hamiltonian the total collective Hamiltonian is first divided into two 
parts, an intrinsic part and a part that only describes the radial motion, 

H = -—R^dnR'^dR + H{R). (17) 

The hamiltonian H{R) is nothing but the sum of the potential v{R) we used in the pertur- 
bative calculation, plus the kinetic term used to calculate the energy difference in the energy 
denominators, 
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We now introduce the instantaneous eigenstates of the intrinsic Hamiltonian. For fixed 
R we diagonahze if, 

H{R)\^^ = E.,{R)\^R- (19) 

The standard Born-Oppenheimer approach is to start at large i?, where the states are pure 
A^A^, and follow the continuous energy curve as R becomes smaller to define the adiabatic 
potential energy. This is basically the procedure we follow, apart from a few additional 
complications encountered if we have more than one A^A^ state in any given channel. This 
is solved by using a model Hamiltonian to determine the unmixed adiabatic eigenstates, as 
will be discussed below. 

Having defined the adiabatic eigenstates, the mathematically correct way to derive the 
adiabatic Hamiltonian is to decompose a general time-dependent solution of the Schrodinger 
equation in adiabatic eigenstates, 

\^>{t)) = Y.aAR,t)\€)^. (20) 

i 

If we substitute this in the time-dependent Schrodinger equation, and multiply from the left 
with ji(i\^ we find 

j 

+ E R{A^j-,dRR'dn\j)^a,{R,t). (21) 

j 

This expression exhibits both the adiabatic Hamiltonian (the first term on the right-hand 
side), and the corrections. If we assume there are no diagonal terms of dn between adiabatic 
eigenstates (since there are no closed circuits in R space, it seems a reasonable assumption to 
ignore Berry's phase), the first correction describes a velocity dependent transition potential. 
The next term gives a non-adiabatic correction to the diagonal adiabatic potential, as well 
as a velocity independent transition potential. 



1. The potential in angular momentum coupled form 

To evaluate (|IT]) it is convenient to introduce states \I1I2LSJT) labeled by the conserved 
(total) isospin T and angular momentum J, as well as the NN channel orbital angular 
momentum L and spin S. Since we wish to evaluate the matrix elements of H{R) in the 
basis \I[l2LSJT) we need to know the matrix elements of the potential v in this basis. After 
performing some standard angular momentum algebra, we find that 

{hhLSJTl V \I[I'2LSJT) 



vi5ss'hL'Si^i'^5_ 



hi'2 
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+|(-ir%5'5LLJ;:;:^}|;:;:^} (/iiiit(-) i/o (/,iiit(-) %) 




t(^^) \i[) (/2|||t(^'^ 1/2) • (22) 



The relevant reduced matrix elements of T are 



(Ar|||T(") \N) =-10, (23) 

(A|||T(^i) I A) =-20, (24) 

(A|||T(") \N) =-8^2. (25) 



The matrix elements of the kinetic part are taken to be very simple, 
{hhLSJT\ K \I[I'^L'S'JT) = 6i,r^5j,j^6LL'Sss' UoO[h + h- 1/2] 



L{L 



2Mj,j,R\ 



(26) 



To account for the effects of centripetal repulsion at small R, we have added the centrifugal 
energy, h^L{L + l)/2Mi^i^R^^ to channel energies. We take Mi^i^ to be the reduced mass in 
the relevant channel, 

where M1/2 = 932 MeV and M3/2 = 1232 MeV. 

For sake of comparison we shall also need the matrix elements of a model nucleon-nucleon 
interaction in the same channel (we have suppressed the spin and isospin projection quantum 
numbers), 



{NNLSJT 


e + V,^a' 


■a^ + Vi 


V,V|(3i?, 


i?,- 


^^^) 


NNL'S'JT) 


= ^SS'^LL' 




-l)i+^6< 


1/2 1/2 s 1 
1/2 1/2 1 * 


+ 








6V3ov;^( 


'_l^L+L'+S+Jll,gg, 


(l2L'\{ 


L S 
S' V 


J 

> < 

2 


'1/2 

1/2 

1 


1/2 s " 
1/2 s' ■. 
1 2 



(28) 



Here we have chosen to use a parametrization identical to the one used in the discussion of 
perturbation theory. 

In order to determine an effective potential from the analysis given above we perform a 
calculation that is standard in the Born-Oppenheimer approximation. We assume that at 
each R we can diagonalize the potential. We then follow a state from large R, where it is 
purely made up out of nucleons, to small R where we have mixing. This gives the usual 
adiabatic potential curves of molecular physics. A complication will be that we may have 
more than one nucleon state for given quantum numbers. We shall see that we can still fit 
the interaction in that case. 
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The model space splits into two parts that are either symmetric or antisymmetric under 
interchange of the particles [^]. We select the antisymmetric part where the NN states 
satisfy the standard selection rule L + S" + T is odd. Of course, in contrast to the lowest 
order perturbative result, we will now find a channel dependent potential. 

The fact that we have a different reduced mass in each channel will also lead to a position 
dependent effective mass. Suppose that cj_^j^{R) gives the percentage of the /1/2 state in 
the diagonalization of the Born-Oppenheimer Hamiltonian. In that case the average radial 
kinetic energy is 

y ^cj jiR)\dRR^dR, (29) 

and we can define a position dependent mass by 

M{R) = 1/ f E %^1 • (30) 

\hI2 ^^^hh 



2. Derivation of the potential 

Let us first consider the case T = 0. A little analysis shows that in this case the two 
channels J'^ = 1~^ and J'^ = 1^ give three two-nucleon states, enough to determine the three 
independent functions in (pS]). (There are no A^A^ states with T = and J = 0.) 

For J'^ = l"*" we have the two possible states \NNL = OS = IJ = IT = 0) and 
|A^A^L = 2S = IJ = IT = 0). The matrix elements of (E8[) in this space take the form 



2v^rt K° + K° 




(31) 



with eigenvalues 

El = V;° + V^° + 2V,°, E2 = K° + K° - 4Vt°. (32) 

These numbers should be fitted to the two lowest eigenvalues of the diagonalization of 
the algebraic potential in the complete nucleon-delta space. If we also look at the state 
|A^A^L = 1S = 0J = 1T = 0), with matrix element of (^) equal to 

E, = K° - 3V;°, (33) 

(again, this is equated to the lowest eigenvalue of the complete diagonalization) we find 

V;° = Ei/2 + E2/A + E3/4, (34) 

V^ = Ei/6 + E2/I2 - E3/4, (35) 

V,' = {E, - E2)/6. (36) 

A similar calculation for T = 1 involves the three channels J'^ = 0"'",0^,1~ (again 
L + S + T =odd). The relevant matrix elements of ( p8|) are 
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E[ = {NNL = OS = OJ = OT = l\V \NNL = OS = OJ = OT = 1) 

= K' - 3K' 

E'^ = {NNL = 1^ = IJ = IT = 1| l^ \NNL = IS = IJ = IT = 1) 
E'^ = {NNL = IS = IJ = OT = 1\V \NNL = IS = IJ = OT = 1) 

= \/' + K'-4v;^ (37) 



Thus 



V,' = E[/A + E'j2 + EyA, 

V^ = -E'JA + E'JQ + E'J12, 

V,^ = {E', - E',)IQ. (38) 

We subtract the centrifugal force from our potentials, since in general this is not taken to 
be included in the potential itself. We had to include it in our kinetic term, since we make 
the Born-Oppenheimer separation between the radial motion and everything else, including 
the orbital motion. 

IV. RESULTS 

As stated before, for each total isospin, T = 0, 1 we write V in the \NNLSJT) space as 

{NNLSJT\V^ + V^a^ ■ a^ + V^a](7]{'iRiRj - 6ij)\NNL'SJT). (39) 

It is the V^ , V^ and V^^ that we study. We have extracted the necessary values of f i, U2, v^ 
from the work of Walhout and Wambach 0. This is an uncertain process for large R 
since the potentials are quite small in this region and are calculated in P] on a lattice by 
first calculating the total interaction energy of the B = 2 system and then subtracting the 
rest energy of two free Skyrmions. We wish to compare the Skyrme based potential with 
a realistic nucleon-nucleon interaction. One cannot relate our result to modern effective 
potentials, such as the Bonn [^] and Paris ||25| potentials, since their central parts contain 
explicit momentum dependent terms. For that reason we compare our potentials to the Reid 
soft core (RSC) potential p6[| . 

In Figs. 1 and 2 we show the central potentials V^ calculated from Eq. (|TT|), first from 
the first term on the right of (p!T| ) only (this is the result of Ref. [^) and then with the 
Born-Oppenheimer method. Note we only show the result for intermediate i?, 1 fm < 
R < 2 fm. To keep the figures simple we separate T = and T = 1. The nucleons only 
result from is independent of T and shows a weak attraction. The perturbation result 
(p!5|) gives considerably more attraction, which is in turn somewhat reduced by the full Born- 
Oppenheimer (BO) diagonalization. This reduction comes mostly from the centrifugal terms 
in the kinetic energy which were not included in the perturbative calculations, and partly 
from the higher order effects included in the diagonalization. The centrifugal terms reduce 
the mixing at small R, and hence lead to less central attraction. All the results, nucleons 
only, perturbation theory and BO diagonalization agree at large R. This reflects the fact 
that state mixing involves multiple pion exchange and hence has intermediate range. Also 
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shown in Figs. 1 and 2 are the T = and T = 1 components of the Reid soft core potentiaL 
They have less central attraction than we find. As we shall see that reflects a rather different 
shuffling of attraction between the central and the spin dependent interactions between the 
RSC and the Skyrme approach. But even at this level the fact that the combination of careful 
solution of the B = 2 sector, state mixing and finite Nq effects gives too much attraction is 
a welcome change from the many previous attempts to derive the A^A^ interaction from the 
Skyrme approach which always found too little central attraction. 

In Figs. 3 and 4 we show the T = and T = 1 spin-dependent potentials. In these 
cases the nucleons-only potential, the perturbative results and the full diagonalization are 
all quite similar, but now the RSC is stronger than the Skyrme result, rather than weaker as 
it was in the central potential. This is the shuffling we referred to above. Figs. 5 and 6 show 
the tensor force for both values of isospin. Except for T = 1 at small R all the Skyrmion 
calculations and even the RSC agree. This reflects the dominance of one pion exchange for 
the tensor potential. 

We see that though the effect of the Born-Oppenheimer corrections are modest for large 
R they become large for intermediate R in the central channel, but remain small in the other 
channels. The reason for this is that the corrections are small for all channels compared to 
the typical size of the potential, but that the lowest order central potential is much smaller 
than this typical size. 

We saw in Figs. 1-4 that the RSC gives less central spin-independent attraction than we 
flnd but a stronger spin-dependent potential. It is therefore instructive to compare the more 
physical S'-channel L = potentials. In Fig. 7 and 8 we show the ^5*1 (T = 0) and the ^5*0 
(T = 1) potentials. Now we see that the full Born-Oppenheimer potential looks quite similar 
to the RSC while the one obtained without mixing does not resemble the RSC. It is in this 
sense that we say that the present results give a credible account of the nucleon-nucleon 
interaction at intermediate range. 

We now turn to a more detailed analysis of the Born-Oppenheimer approximation, the 
adiabatic approximation and non-adiabatic corrections. We analyze a representative case, 
the J = 1~, T = 1 channel, since it mixes A^A^, A^A and A A states. The relevant states 
in the diagonalization are |A^A^, L = 1, S" = 1), |A^A, L = 1, S" = 1), |A^A, L = 1, S" = 2), 
I AA, L = 1, S* = 1) and | AA, L = 1, S = 3). In Fig. 9 we plot the energy of the adiabatic 
eigenstates, which remain distinct for all separations, except for a narrow avoided crossing 
near 1.2 fm. For this channel all the energies are repulsive. In Fig. 10 we show the percent- 
age of each component in the lowest adiabatic eigenstate. It is this state that was used in 
the determination of the A^A^ potential. One can see that the A^A^ component is dominant, 
and we have only small admixtures of the other components. For small R it is the third 
adiabatic eigenstate that admixes most strongly. 

In Fig. 11 we show the effective mass of the lowest adiabatic eigenstate as deflned by 
Eq. (pop. Since the only effect on the masses we include is the centrifugal term, the mass 
increases with decreasing R. Preliminary calculations of the mass parameters for the classical 
B = 2 system give mass parameters that decrease with decreasing R, at least in the attractive 
channel |2^. Thus the results of Fig. 11 should be taken only to remind us that the Skyrme 



approach leads to a picture of nucleons that distort and change their mass as they interact. 

Corrections to the adiabatic approximation lead to transition potentials among the adi- 

abats. The relevant matrix elements can be calculated exactly using standard algebraic 
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techniques: 

R {t\ Or \j)j, = R{t\ (dnHiR)) \j)^ /{E,{R) - E,{R)). (40) 

We further assume that the (imaginary) diagonal matrix element is zero, which corresponds 
to a consistent real choice of |j)^. The second derivatives can now be evaluated using a 
sum over intermediate states. As argued in the previous section we have two terms, a local 
transition potential and a velocity dependent potential. We continue to study the J=l~, 
T = 1 channel and consider the transition from the A^A^ adiabat to each of the five others. 
The local potentials for these transitions are shown in Fig. 12 as a function of R. They are 
on average much smaller than the difference between the relevant adiabatic eigenvalues, as 
they should be if the adiabatic approximation works. Near R = 1.2 fm there is a spectacular 
excursion of the velocity dependent potential in Fig. 12. This can be traced back to the 
narrowly avoided level crossing in Fig. 9. This picture would suggest that at high A^A^ energy 
there should be the onset of strong A production. The theory is still too crude to make 
a detailed prediction, but this aspect of A^A^ dynamics suggested by the Skyrme approach 
deserves further investigation. The velocity dependent potential are a little more difficult to 
interpret. We have chosen to plot, in Fig. 13, the potential such that we have to multiply it 
by —i times the scaled radial velocity —ip = v/c = -^Qr. As can be seen this potential is 
very small for non-relativistic velocities. 

V. SUMMARY AND CONCLUSIONS 

Here we have shown that the Skyrme model can give strong mid-range nucleon-nucleon 
attraction that is in qualitative agreement with phenomenological potentials. This is in sharp 
contrast to early results based on the product ansatz. Two components play a key role in 
obtaining the attraction. First one must pay attention to the non-linear nature of the Skyrme 
Lagrangian, as shown in Refs. P,|14| , but second one must also include configuration mixing 



at intermediate distances as we have stressed here. This mixing is not easily formulated 
in the large A'c limit, but is easily included for A'c = 3, and hence brings with it some 
other finite A'c effects. The combined effect of the careful treatment of the non-linear 
equations and the configuration mixing is to give substantial central mid-range attraction 
for the NN system that is in qualitative agreement with the data. To go from this work 
to a theory that can be confronted with experiment in detail is a difficult challenge. There 
are R dependent corrections to inertial parameters to include, there are dynamical quantum 
corrections, there are non-adiabatic effects that are particularly important at small i?, and 
there are other mesons to include in the Skyrme Lagrangian. All these are under study. The 
success so far encourages us to proceed and to suggest that the results obtained so far will 
be robust under these refinements. 

The success of the Born-Oppenheimer approach to the nucleon-nucleon interaction in the 
Skyrme model raises interesting questions about what "nucleons" are when they are interact- 
ing. We find that the state that is two nucleons at large separation becomes, adiabatically, a 
state with delta components as the nucleons approach. For the Skyrmion-Skyrmion system 
itself this distortion and deformation under interaction is even more dramatic. What then 
should we say about the interacting nucleons? That there is a big delta component in their 
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wave function? That they are distorted by the interaction? In some sense this is true, but 
in the spirit of the Born-Oppenheimer approximation, it is not. The adiabatic potential 
derived from the Born-Oppenheimer prescription involves state mixing, in fact much of its 
attraction comes from that mixing, but that adiabatic potential is to be used as an effective 
potential between nucleons. Think of the van der Waals potential between hydrogen atoms. 
It comes from virtual excitation of the atoms, principally to the first L = 1 state, but the 
potential is to be used between atoms in their ground state. One usually does not talk about 
the percent of atomic L = 1 state in the hydrogen molecule. Following that example, we 
should not talk about the percent of delta in the state of two interacting nucleons. On the 
other hand, one could take apart the Born-Oppenheimer calculation and quote such a per- 
cent. It is really a matter of what basis one uses. Many discussion of delta's in the nucleus 
are unclear about this ambiguity. Our calculation helps to emphasize the perspective de- 
pendence of these discussions. We plan to return to this topic both for the nucleon-nucleon 
scattering system and for the deuteron in a subsequent paper, as well as to the question of 
the interaction effects on the nucleon mass parameter. We have also seen that corrections 
to the adiabatic approximation make it possible to compute transition potentials. These 
can be used to calculate Delta production in A^A^ collisions. In this paper we have only 
scratched the surface of this interesting problem. 

In summary we have shown that including the effects of channel coupling (A mixing) 
in the projection of the nucleon-nucleon interaction from recent studies of the potential 
in baryon number two Skyrmion systems, substantially increases the strength of the mid- 
range central attraction bringing it into qualitative agreement with experiment. To obtain 
this result we introduce some quantum corrections (equivalently corrections for the finite 
number of colors). Many more sophisticated quantum corrections remain to be made. But 
our results show that a non-perturbative approach to the problem of obtaining the "static" 
nucleon-nucleon interaction from QCD based on the Skyrme approach has great promise. 
We are therefore encouraged to follow up the first success reported here. 
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APPENDIX A: DETAILS OF THE PERTURBATIVE EXPANSION 

We exhibit details of the calculation for the first matrix element, for the part of H 
proportional to W. We again make use of an equation from |jl5|, Eq. (30). 



{NN\ W \NA) {NA\ W \NN) 

= {N\ T- \N) {N\ T- \N) {N\ T^, | A) (A| T,^- |Ar) 

= {N\ T^, \N) {N\ T^^ \N) {{N\ T^%^^ \N) - {N\ T^, \N) {N\ T^ \N)) 

= {NN\ -fr^alria,PUlQU|S^, ' ^^^i^lWv, ' fW^')} 1^^) 
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= {NN\ N^^^{2 + fa^- a^){2 + fr^- r^) \NN) . (Al) 

Here we have used the finite- A'"c correction terms Pat, defined before, and Qn = 



1 — 1/A^c)(l + 5/A^c)- In this manner we find that 

{NN\ H{R) \N^) (A^AI H{R) \NN) 



3' "' ;i9V'^2 "T ^^3; "T 27V'^2 ~ "^s)^ ■ ^ 



= {NN\ P^g^(2 + It' . t'){\{vI + 2vl) + ^^{vl 

1y(^3 - V2Vz){'ia' ■ Ra^ ■ R - a' ■ a^)} \NN) , (A2) 

and for an intermediate state with two A's, 

{NN\ H{R) |AA) (AA| H{R) \NN) 
= {NN\ Q%{1 - Y . t'){\[vI + 2vl) - ^^[vl - vl)a'- a' - 



m 



V2V3){3a'- Ra^-R- a'-a^)} \NN) . (A3) 



When we use the fact that the energy needed for excitation of a single A, which we denote 
by 6, is half that for excitation of two A's, we find 

KLI = -^{[(Ip'n + IQI) + {4fPN - ^Qiy- r\vl + 2vl) + 
{a'- a^)[(^P^ - ^Ql) + (^P^ + ^,Q%y- r']{vl - vl) 
iSa'. Ra'. R-a'. a^)[(^P^ - ^Q^) + (^P^ + ^,Q%y- r^Kv^ - v,Vs)}. 

(A4) 

Of course this can be rewritten by introducing spin and isospin projection operators, 

pr_ l-n-T-2 

PI = ^^, (A5) 

ri.r' = -3P([ + P;. 

We thus get the form (subscripts now denote isospin projections) 

vSL = -^miPo + im + ifQDPiKvi + 2vi) + 

y-<^')[-T,Q"NPo + {f,PN - ^,Ql)PlM - vl) 

{3a'- Ra'- R-a'- a')[-±Q%P^ + (ifp^ - ^Q^)P^-](^;2 _ y.y,)}. 

(A6) 
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FIGURES 

FIG. 1. The central potential V^ as a function of R in the region 1 to 2 fm for the T = 
channel. The solid line gives the nucleons only result of |^. The longer dashed line is the result of 
the state mixing in perturbation theory and the shorter dashed line of the full Born-Oppenheimer 
diagonalization. The dotted line is the Reid soft core potential in this channel. 



FIG. 2. Same as Figure 1 but for T = 1. 

FIG. 3. The spin dependent potential V^ as a function of R in the region 1 to 2 fm for the T = 
channel. The solid line gives the nucleons only result of |^. The longer dashed line is the result of 
the state mixing in perturbation theory and the shorter dashed line of the full Born-Oppenheimer 
diagonalization. The dotted line is the Reid soft core potential in this channel. 

FIG. 4. The same as Figure 3 but for the T = 1 spin-dependent potential. 

FIG. 5. The tensor potential V^ as a function of R in the region 1 to 2 fm for the T = 
channel. The solid line gives the nucleons only result of |^. The longer dashed line is the result of 
the state mixing in perturbation theory and the shorter dashed line of the full Born-Oppenheimer 
diagonalization. The dotted line is the Reid soft core potential in this channel. 

FIG. 6. The same as Figure 5 but for the T = 1 tensor potential. 

FIG. 7. The potential in the ^Si (T = 0) channel as a function of R. The solid line gives the 
nucleons only result of Q. The dashed line is the result of state mixing in a full Born-Oppenheimer 
diagonalization, and the dotted line is the Reid soft core potential in this channel. 

FIG. 8. Same as Figure 7 but for the ^Sq (T = 1) channel. 

FIG. 9. The Born-Oppenheimer energy eigenvalues coming from diagonalization of the adia- 
batic energy in the J = 1^, T = 1 channel as a function of the adiabatic variable R. The lowest 
energy state goes over to the A^A^ state for large R. The next three states go over into the A^A 
states, and the upper two into the AA states. 

FIG. 10. The percentage admixture of configurations in the lowest {N — N ) adiabatic state 
of the J = 1~ T = 1 channel as a function of R. The solid line is the A'^A^ component, and the 
seecond largest component, given by the short-dashed line is the state | A^A, L = 1,5 = 2). 

FIG. 11. The effective mass of the "nucleon" in the lowest J = 1^ T = 1 channel as a function 
of R. This is taken to be twice the reduced mass. 
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FIG. 12. The velocity independent transition potential from the lowest (NN) adiabatic state 
of the J = 1~ T = 1 channel as a function of R for transitions to each of the five other channels. 
The numbers in the legend refer to the order of the states in Fig. 9, where one is the lowest energy 
adiabatic state. 

FIG. 13. Same as Figure 12 but for the velocity dependent transition potential. Each term is 
understood to be multiplied by —i(3 where f3 = v/c. 
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